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The Drude weight for the one-dimensional Hubbard model is investigated at finite temperatures 
by using the Bethe ansatz solution. Evaluating finite-size corrections to the thermodynamic Bethe 
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■ ansatz equations, we obtain the formula for the Drude weight as the response of the system to an 

' external gauge potential. We perform low-temperature expansions of the Drude weight in the case 

^v^j , of half-filling as well as away from half-filling, which clearly distinguish the Mott-insulating state 
from the metallic state. 
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I. INTRODUCTION 



The Mott-Hubbard metal-insulator transition (MIT) is one of the long-standing important issues in strongly corre- 
' lated electron systems. The one-dimensional (ID) Hubbard model is ^fundamental model which describes the MIT. 
This model is exactly solvable in terms of the Bethe ansatz methodlil. Various thermodynamic quantities such, as 
the specific heat and the spin- and charge-susceptibilities which characterize the MIT have been obtained exactlyBa. 
However, the study on transport properties such as the Drude weight based upon the Bethe ansatz method was 
restricted to the zero tcnu)|Urature caseETLl. The Drude weight at finite temperatures-iiKas investigated only by using 
numerical methods so faiEra. Since the Drude weight is a direct probe for the MITLL9eI, it is desirable to obtain the 
^ ■ exact results on it at finite temperatures. This is the main purpose of this paper. In order to obtain the Drude weight, 
it is necessary to calculate the finite-size corrections to the energy spectrum, so that we should evaluate the finite-size 
corrections to the thermodynamic Bethe ansatz equations. For this purpose, we generalize standard methods for 
finite-size corrections based on the Euler-Maclaurin formula to the case of finite temperatures, and obtain the leading 
, temperature dependence of the Drude weight at low temperatures. 

■ It is easily seen from the Kubo formula that the Drude weight in translationally invariant systems is independent 
^\ \ of temperature. Then the temperature dependence of the Drude weight for the Hubbard model comes from the 
"j***; ■ interactions such as the Umklapp scattering which break translational symmetry. This implies that its temperature 

dependence carries the information how the electron correlations which control the MIT develop in low temperatures. 
In the case of half-filling, we indeed show how the system approaches the Mott insulator as the temperature is 
decreased. 

The organization of this paper is as follows. In Sec. II, the formulation for the finite-size corrections to thermody- 
namic Bethe ansatz equations is given. The procedure is not restricted to the Hubbard model, but also applicable to 
any other solvable models. In Sec. Ill, we obtain the general expression for the Drude weight at finite temperatures. 

■ The low-temperature expansion for the Drude weight is derived in Sec. IV both in the case of half-filling as well as 
away from half-filling. The leading temperature dependence is evaluated at low temperatures. Brief summary is given 
in Sec. V. 



II. FINITE-SIZE CORRECTIONS TO THE THERMODYNAMIC BETHE ANSATZ EQUATIONS 

In this section we consider finite-size effects on the Bethe ansatz solutions of the ID Hubbard model at finite 
temperatures. Finite-size effects at zero— temperature have been studied by many authors in connection with the 
application of conformal field theorytlTo. We generalize their method to the case of finite temperatures. The 
hamiltonian of the ID Hubbard model reads, 

H = -^2 4^+1 + h.c. + c{i c fi c ii c ^- w 
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It is necessary to introduce the Aharonov-Bohm (AB) flux $ to formulate the Drude weight as the response to an 
external gauge potential. Alternatively, the.-effect of the AB flux is incorporated into the twisted boundary condition 
for n thje wave function, ty(x + L) = e J *5 , (x)t-3. The Bethe ansatz equations in the presence of the AB flux are given 
byllB, 

e ikjL = e i# TT sin k 3 -A a + iu 
*■ sin ki — A„ — iu ' 

a=l 

N . . , , . M 



nl\ a — Sill Kj -f IU T-r ll a — l\p -f Zi'U 

A a - sin fc, ■ - iu ~ *■ A a - Ag - 2iu ' 

j=i a 0=1 a p 



Here kj and A Q are, respectively, the rapidities for the charge and spin degrees of freedom, and we have introduced 
u = U/A. N and M are the total number of electrons and down spins. 

Thermodynamic Bethe ansatz solution for the ID Hubbard model was obtained by Takahashi many years ago with 
the use of so-called string hypothesise!. The validity of the string hypothesis is justified only for the thermodynamic 
limit L — > oo. In order to calculate the Drude weight, we should evaluate the energy for a finite-size system, because 
the effect of $ vanishes in the thermodynamic limit. Thus one may worry about whether the string hypothesis can 
be applied to the calculation of the Drude weight. However, by recalling the following fact we can still adopt the 
string hypothesis for our purpose: The corrections to the string hypothesis for a finite-size system is estimated as 
~ 0(e~ cL ), where c is a constant which depends on the temperature. On the other hand, the dependence of the 
energy on the AB flux appears in the order of 1/L 2 . Thus the correction to the string hypothesis for the finite-size 
system is much smaller than the finite-size corrections to the energy spectrum which we need for the calculation of the 
Drude weight. This situation is analogous to that for the Kondo modeO or the impurity Anderson mdeOE3, where 
the local electron correlations, which are given by the 1/L-corrections to bulk quantities, are correctly evaluated based 
upon the string hypothesis. 

Using the string hypothesis, we can thus write down the thermodynamic Bethe ansatz equations. After taking the 
logarithm of the above equations, we end up with 

oo M„ 



n— 1 a— 1 

ML 



n—1 a—1 v / 



L(sin" 1 (A"^ + inn) + sin _1 (A'™ - inu)) = 2nJ' n a + 2n$ 

N-2M' / • , a/"\ 00 /A'™ A' m 

sin ki — A' \ v-^ v-^ „ / A - A 



e •(=^i i )+i;i;M ;i ? i ). 



N-2M' 



sinfc,-A«\ . „ . A _^ /AS -A 



j=l v 7 m=l v 7 

with 9{x) = 2 tan -1 x and 

-^-)+29( ] ^—)+... + 28( ?—S) + o(-E—) n ^ m 

Q nm {x) = { %\ m V V"x m| > J fx\ Kn + m - 2j Kn + mJ (7) 
-2) + --- + 2d {—2) +0 {^) 

Here kj is the rapidity for charge excitations which are not in bound states, A™ is that for spin excitations, and A'™ is 
that for bound states. Ij, J™, and J' a are the corresponding quantum numbers which specify the above excitations, 
respectively. M n is the number of n-strings for spin excitations. 2M' is the total number of electrons which make 
bound states. 
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In order to calculate the finitfcsize corrections to the energy spectrum, we expand the rapidities in terms of 1/L 
following Berkovich and Murthyo, 

k s =k? + (± + ^ + 0(1/1% 
A «=Ar + ^ + ff + 0(1/L% 

A /n =A/ „oo + /Hn + /^ + 3 
L L A 

The lowest-order contributions in 1/L give the conventional thermodynamic Bethe ansatz equations which readi, 

<-\^rtuw fu\ 1 _lV^ f°° dA nucosk(a n (A)+a' n (A)) 

(1 + C(*))P(*) = ^ + L V (n U )'-Kainfc-A)' ' (9) 



where 



7 7 „(A) ( T„(A) + ]T A™ * <Tm(A) = f — -, (10) 
^ Joo k {nu) 2 + (sink - A) 2 

v' n (A)a' n (A) + J2 A nm * </ m (A) = -Re 1 , = - /'°° - . A , 2 , (H) 

^ 77 v/1 - (A - mu) 2 Joo tt {nu) 2 + (smk - A) 2 



lnC(fc) = T h'-oo ~ (»«)' + (sin ib-A)» ' (12) 



, . .... /" r dk nwcosfcln(l + C 2np H . , 

In 1 + ri n (A)) + / v . s v ;; = — £L_ + ^ A nm * In 1 + r,^ A , 



m— 1 



,'" r dfc nucosfcln(l + C 1 (fc)) 

i °( i+ ^(A))+ 1 - w+( ; infc _ A y j 



4rWl - (A - inuf - 2nA ^ „ , , ,-1,,^ 
^ ^= * + J] A nro * ln(l + 77 m (A)), 



(13) 



(14) 



d f°° dx' x — x 
A nm * 4>(x) = 5 nm 4>(x) + — / — 6„ m ( )4>(x'). 

p(k), cr n (A), and cr'„(A) are the distribution functions for the rapidities, kj, A™, and A'™, respectively, and ((h) = p h /p, 

r] n (A) = cr^/(T n , r)' n (A) = u' n Ja' n with the distribution functions for holes p h , a^, and a' n . We have introduced an 
external magnetic field H and a chemical potential A 

Using the Euler-Maclaurin formula and eqs.(^)-(^|) and (§), we obtain the 1/L- and l/L 2 -corrections to the Bethe 



ansatz equations which determine /i,2) <?i, 2n an d ^l,2n- Ta 



dng a continuum limit, we consequently have, 



(l + C(*)M*)/i(*) = ^- + E / — 7 i2-m rxa (gi»(AK(A) + fe lw (Ay TO (A)), (15) 

Z7T n=1 J-oo 77 (.™J -h^Sin/C — IV) 

a n (A)Tin(A)g ln (A)+jrA nm *a m (A)g ln (A)=[ — - ^ cosfc — j 1 ( A; ) p ( fc ) j ( 16 ) 

J—tt 7T ^77,7/ j -f- [SID. /£ — A) 
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a' n {A)r,' n {A)h ln (A) 



OO 

E 



K a i / 1 \ i / \ \ n * f dk 
. A nm * a m (A)/ii„(A) = / 

1 7T J-7T 7T 



! COS fc 



(um) 2 + (sin fc — A) : 



rh(k)p(k), 



(17) 



oo „ c 

(l + C(*0)p(A)/a(*) = £ / 
— i •/— < 



dA 



7T (to) 2 + (sinfc — A) : 



r (g 2 n(A)o- n (A) + ft 2n (A)c/„(A) 



-^[(i + C(fc))p(fc)A 2 (fc)] 



E 



dA to (sin k — A) 
7r ((to) 2 + (sinfc - A) 2 ) 2 



( ff 2 n (A)a„(A) + / l 2 „(A)a' n (A)), 



(18) 



(1 + ^(A))a n (A) 52 „(A) = + ??„(A))(T n (A)g'f n (A)] 



x dfc 



1 cos fc 



_ w tt (™) 2 + (sinfc-A) 2/2(fc)p(fc) 



*• dfc 



1 sin fc 



/i 2 (fc)p(fc) 



w 7r (to) 2 + (sin k — A) 2 



y. /■°° dA' e , . A- A' <b„,(A>„(A') 



m—1 



lim 



° dA'_,A-A'. 5 2 m (A> m (A') 



1 



Aq^oo 487TW 



E 



u 2u 2 
e'((A-Ao)/«) 



w dfc to (sin fc — A) cos 2 fc 
. w 7r ((to) 2 + (sinfc - A) 2 ) 

2A' 
2^ 



e'((A + Ao)/«) 



2 f?(k)p(k) 



(1 + 77 m (A ))cr m (A) (1 + ?7 m (-Ao))cr m (-A) 



(19) 



(l + rj' n (A))a' n (A)h 2n (A) = ^[(l + rj' n (A))a' n (A)hl n (A)] 



71 dk 



nu cos fc 



_ n 7T (TO) 2 + (sinfc-A) 2/2(fc)p(fc) " ./_„ 7T ((nu) 2 + (sinfc - A) 2 ) 2 



T dfc to (sin fc — A) cos 2 fc 



A 2 (fc)p(fc) 



T dfc 



/ 2 (fcHfc) + yW°° ^ 9 , ( ^ ) maw 



7r (to) 2 + (sinfc - A) 2 2 

dA' „ ,A-A'. h? m (A') CT ' m (A') 
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m—1 



/If-' 



1 



lim 

Ao^oo 487TU 



M 2w 2 
0'((A-Ao)/«) 



e'((A + Ao)/«) 



(1 + n' m (A ))a' m (A) (1 + ,/ m (-A ))a' m (-A)) 



(20) 



where 0'(x) and Q"(x) are, respectively, the first and second derivative of 0(x). This gives our starting equations 
for the following discussions. Using these equations, we shall investigate finite-size effects to obtain the Drude weight 
at finite temperatures. 



III. DRUDE WEIGHT AT FINITE TEMPERATURES 



Here we derive the expression for the Drude weight at finite temperatures using the formulation in the previous 
section. The Drude weight at finite temperatures is given by the second derivative of the energy spectrum with respect 
to the AB flux dl2, 



D= L/d 2 E ri 



2 \ d<f> 2 



(21) 



6=0 



Here (• • •) is the thermal average for a canonical ensemble. Note that the above Drude weight D is different from the 
second derivative of the free energy, which denotes a Meissner fractionEd. 
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We now wish to evaluate D from the finite-size corrections to the energy. To this end, we first write down the total 
energy, 

„ N—2M' oo 

- = - (cos k 3 + fioH + + ^V 1 - ( A '2 - i™) 2 

j — 1 n—1 a. 

oo oo 

n—1 n—1 

We then expand the energy in powers of 1/L using eq.(|8|), 

E Ei Ed 

T -E + T + -. (23) 



The first-order correction term E\ is given by 



/TP OO /• OO /A ■ \ 

dfcsinfc/ 1 (fc)p(fc) + ^ / dA4Rc— F =L=S=/ii n (A)cr'„(A). (24) 
-7T n= l"'-00 V 1 



(A — inu) 2 



Differentiating eqs.(|12|)-(|l4j) with respect to rapidities, and substituting them into eq.(24), we easily find that E\ = 0. 
We now consider the second-order term which is, 



E 2 = 2 [ dk sin kf 2 (k)p(k) + 2 f dk cos fe MjMjjj 
f] /°° dA4Re ~ {A ~ mU) h 2n (A)a' n (A) 



n=l ' 



^ /.OO 

„ =1 J-oo 



-1 



<(A)a'„(A) 



(1 — (A — inu) 2 ) 3 / 2 

Using eqs.(^2])-(|20|), we can rewrite this expression as, 

F =- r dk ^ fc )A 2 (fc) f d((k) \ 2 rf> A 00 <7»(A) gl 2 „(A) / dffa(A) 
2 2./_ 7r CW(l + C(fc))l * y 2 il^oo Vn(A)(l + Vn(A))\ dA 

o.Z!v r , A ^(A)feL(A) ( dn' n {K) y 
2 ,tt^oo u' n (A)(l+r,' n (A))V rfA ; ' 



(25) 



(26) 



Note that the dependence on the AB flux appears only through /i(fc) ! gi„(A), and hi n (A). We can obtain the Drude 
weight by differentiating E 2 twice with respect to <£>. We see from eqs.(Eq)-(0) that the equations for dfi/d&, dg\ n /d<& 
and dhi n /d§ do not depend on Thus, we have 

d 2 fi = d 2 g ln = d 2 h ln = 

rf$ 2 rf$ 2 1 ' 

Then the Drude weight is given by 
ld 2 E 2 

dh\ 2 d ( -1 \ 1 dn(k) 



D 



2 d<$> 2 



dfc^(i + c(fc))p(fc)- d$i dfc y 1 + eK «/T; ( i + c(fc))/9 ( fc) dfc 

■ ' ' -1 \ 1 de n {A) 



n—1 v y v 

I|£,a{ (1 + ,/(A)K(A)^} 2 A(- 



+ e s„(A)/Ty (i + J7 (A))<7(A) dA 

-1 \ 1 de'„(A) 



+ e e '-( A )/ T y (1 + 77'(A))<t'(A) dA 



(28) 
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Here we have used the conventional notations, n(k) = T\n((k), e n (A) = Tln7y„(A), and e' n (A) = Thi7/ n (A). In 
order to simplify the expression, it is convenient to define the following quantities, 

£ c (fc) = 27r(l+C(k))p(A0|jp (29) 
e„(A) = 27r(l + »7(A)MA)^, (30) 
6„(A)^2 7 r(l + r ? '(A)K(A)^, (31) 

^ )e (i + C(1mJ ' (32) 

^ S (l + ^(A)^ A) | (33) 

^(^(TTT^T- M 

These quantities have simple physical meanings: £ c , £„„ and £& n correspond to the dressed charges generalized to finite 
temperature. v c , v sn , and vt, n are the velocities for charge excitations, spin excitations, and bound states, respectively. 
Consequently, we end up with the simple formula for the Drude weight expressed in terms of these quantities, 



dA d ( -1 



At: dk V l + e e'»(A)/T 



^ n (A) Vbn (A). (35) 



Note that at finite temperatures not only the charge degrees of freedom but also the spin degrees of freedom contribute 
to the Drude weight. The above formula is the one of our main results in this paper. __ 

To conclude this section, we check that the above formula reproduces the known results&El by taking zero temperature 
limit T — > 0. Since e„(A) > for n — 2,3, ... and e'„(A) > for n — 1,2, 3..., the contributions from spin excitations 
with n > 1 and bound states to the Drude weight vanish for T — > 0. Moreover from eqs. dH) and @ we have 



^ s i(±B) = for T — > where ±B is zeros of £i(A). Thus only the contribution from the charge degrees of freedom 
to the Drude weight survives, 

D = JdK J2 *(«) g<y> = *fc (36) 
zeros of « 

Here K c = £^(fco)/2, and v c = v c (ko) with K(±fco) = 0. Then we reproduce the well-known result for zero temperature. 

IV. LOW TEMPERATURE EXPANSION 
A. Case of half-filling 

In this section, we explicitly derive the temperature dependence of the Drude weight in the case of half-filling at 
low temperatures. In this case, the system is in the Mott insulating state with the charge excitation gap. Thus we 
immediately find D — at zero temperature. However, at finite temperatures it can have finite values as we will see 
momentarily. We consider the case that 2u — A > 0, < 2(Vl + u 2 — u), and 2 — [iqH — A < 0, i.e. the charge 
excitation is gapful, whereas the spin excitation is gapless. In order to obtain the temperature dependence of the 
Drude weight, we need to take a thermal average over a canonical ensemble and consider the temperature dependence 
of the chemical potential A. However, in the presence of the Mott-Hubbard gap, the temperature dependence of 
A appears only through that of the Mott-Hubbard gap, which gives a subleading contribution to the temperature 
dependence of the Drude weight, as is shown below. Thus we can safely ignore the temperature dependence of A. 

Following Takahashi's method, we perform low-temperature expansionsa. As a result, we find that at low tempera- 
tures the Drude weight is mainly controlled by the contributions from the charge excitation, the spin excitation with 
n = 1 and the bound-state excitation with n — 1. At low temperatures eq. (p"2|) , ( p^ ) for n = 1, and (14) for n = 1 are 
rewritten as, 



G 



K (k) = -2cosk-» HA + ( B — „ ..^ itAl + dP, (37) 

J — Q 7T li t I Sill rC — iV J 



x(A) = 2 t , H - 4(ReVl - (A - m) 2 - u) - ^ ^ ^ + ^ _ A , )a £i(A') + C 2 T\ (38) 

e'i(A) = 4u - 2A - C 3 Th K ^/ T 9 " + C 4 Tie-( 4tl - 2 ^/ T , (39) 

where 7 = 2 for H ^ and 7 = 3/2 for H = 0. B is defined by the condition £i(±_B) = 0. In the absence of magnetic 
fields, H — 0, 5 — > +00. £i(A) is obtained from eq.(3£j) by using the Wiener-Hopf method. Substituting the solution 
into eqs.(|37|) and (|39|), we obtain n(k) and e'i(A). The generalized dressed charges 6> 61 and 61 are now determined 
by the derivative of eqs.(p^|), (|l6|), and ( |i~7| ) with respect to which are given in the low temperature limit, 

(l + e"W/ r )£ (fc) = l+ / — (61(A) + 61(A)), (40) 

/- Fi (avt\* rffc ucosfc , ,,. /" s dA' 2it ., ... . 

(1 + e-^/ T )6r(A) = j ^ - u2 + {s - mk _ A)2 Uk) - J b - 47i2 + (A _ A02 6i(A)- (41) 

^vawtn , .> C dk ucosk , , /" s <iA' 2u , . . 

(1 + e )6i(A) = 2+1 , —6 fc - / 2 ■ ,. t— 7761(A). 42 

7-7T I" u 2 + (smfc-A) 2 7-s 71" 4w 2 + (A - A') 2 

Solving these equations we have 6(&) — 1 an d 

6l(A) ~ Dl ^f e (2-^H-A+A(B))/T^ (43) 

61(A) ~ D 2 e- e '^ T , (44) 

where 

The distribution functions for rapidities, p(fc) and ci(A) are given by those for zero temperature. Also, er'i(A) is 
estimated as 

a'r(A) ~ C 'V^ e ( 2 -Mff-A+A(B))/T e ^' l( A)/T (46) 
Using eqs.(|35|) and (|§7|)-(|46|), we finally end up with 



D = V T e -A MH /T + Qr Te 2{1-n H-A+A{B))/T\ + e ~(4u-2A) /T e (2-f, H-A+A(B))/T^ ^ 

where Amh = — 2 + {.iqH + A — A(B) is nothing but the Mott-Hubbard gap, and 

. 1 f B dA ua (A) 

(49) 



<*>(A) = r 

J —7T 



7r(A— sin k) 



8ttu cosh 2u 



Here the first term of eq. (|47|) , which is most dominant, comes from the charge degrees of freedom, whereas the second 
and third terms are the contributions from spin degrees of freedom and bound states, respectively. As seen from the 
above equations, the Drude weight vanishes exponentially at half-filling at low temperatures, reflecting the presence 
of the Mott-Hubbard gap. 
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B. Case away from half- filling 



We next consider the case away from half-filling in the absence of magnetic fields, i.e. k(tt) > 0. We first estimate 
the contribution from the charge degrees of freedom. Since we are concerned with a canonical ensemble, we must take 
into account the temperature dependence of the chemical potential A or fco in the case that the number of electrons 



N 
T 



dkp(k) + 2n / dAcr\ 

-* n=l 



(A), 



(50) 



is fixed. At low temperatures eq.(|5C|) is approximated by 

Po(k) 



N 

T 





dk 


J — 7r 




fco 




7T 





1 



oK(k)/T 

A _ i sin fco - 
— tan 

, ir U 



A 



CTl (A) 



ir 2 T 2 8 fpo(k) 



3 dK\K'(k) 



A: = /.'o 



where po(k) is the distribution function for fc at T = 0. Thus we have the temperature dependence of fco, 

- n 2 T 2 8 fp (ky 



Sk 



6/Oo (kn) dK\K'(k) 



(51) 



(52) 



Here fco is fco at T = 0. Then using eq.(|35|), we obtain the contribution from the charge degrees of freedom to the 
Drudc weight at low temperatures, 



K c v c 

7T 

K c v c 

7T 



2tt 12 8k 2 
CT 2 , 



k—kr, 



ttT 2 d 2 
12 8k 2 



(e c (k)v c (k)) 



k—kn 



where the coefficient of the quadratic term is 

a 2 



c 



TV 

12 



8k 2 



8 ( p (k) 



k—kn 



po(fco) dn\K'(k) 



k—k< 



(e c (ko)vc(k )) 



(53) 



(54) 



Here £ c (k) and v c (k) are calculated at T — 0. 

In a similar manner, we can evaluate the contributions from the spin degrees of freedom as well as from the bound 
states. As a consequence we find that they are subdominant compared to the contribution from the charge degrees of 
freedom at low temperatures. Therefore the leading term of the Drude weight is given by eq.(|5^). We expect that the 
expression (|53|) and its temperature dependence ~ -D(O) + CT 2 may be general for all integral models with massless 
excitations. 



V. SUMMARY 



By using the Bethe ansatz solution, we have obtained the formula for the Drude weight of the Hubbard model at 
finite temperatures. The present general formulation is not restricted to the Hubbard model, but also applicable to 
any other integrable models. We have then performed low-temperature expansions both in the case of half- filling as 
well as away from half-filling. In the case of half-filling, the Drude weight decreases exponentially, as the temperature 
is lowered, reflecting the presence of the Mott-Hubbard gap. In the case away from half-filling, it behaves like 
~ D(0) + CT 2 , with the coefficient C expressed in terms of the velocity for charge excitations, the dressed charge 
and their derivatives. Although the essential properties of the Drude weight can be seen through the present low- 
temperature expansion, it is interesting to obtain its full-temperature dependence by solving the integral equations 
numerically, which should be done in the future study. 
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